Coherent Quantum Tunneling of Bose-Einstein Condensates in Optical Lattices. 
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We study inter- well coherent quantum tunneling of Bose-Einstein condensates (EEC's) in optical 
lattices. Irrespective of whether the BEC is a superfluid in the whole lattice, it may or may not 
be a superfluid within each well of the lattice depending on whether the ratio of the well kinetic 
energy to the interaction energy is small or large. When the bosons in each well are superfluid, we 
have Josephson tunneling among the wells driven by the phase difference between the superfluids 
in the different wells. For typical lattice parameters the Josephson tunneling rate will be somewhat 
enhanced from the naive estimate of (tunneling energy scale)/fi. Whether it is enhanced or not 
will depend on lattice parameters, but it can never be zero. When the bosons in each well are not 
superfluid, there is coherent quantum tunneling of essentially free bosons between the wells. This 
tunneling is proportional to the boson number difference between the wells. 
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In the present note we will study Bose-Einstein conden- 
sate (BEC) tunneling in an optical lattice. It has been 
verified experimentally in one dimensional optical lattices 
that there is coherent quantum tunneling in BEC's0. 
We will try to take an in-depth look at tunneling in op- 
tical lattices in 3 dimensions. We take the well tunneling 
amplitude, or hopping, to be t. In order for the tunneling 
to be relevant, this energy scale needs to be larger than 
the kinetic energy scale for the whole lattice: 
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where m is the mass of the bosons and fl is the volume 
of the whole lattice. If the kinetic energy of the lattice 
overwhelms the hopping, the tight-binding approxima- 
tion with bosons hopping from lattice site to lattice site 
with tunneling t is no longer valid. 

On dimensional grounds one expects the tunneling rate 
to be proportional to t/h. On the other hand, the pro- 
portionality factor should be a universal dimensionless 
function of the number of lattice sites Nl, the number 
of bosons N and/or the number of bosons in the conden- 
sate No, and the different energy scales in the problem. 
Besides t and fi^ /2mil?^^ , there is the on-site repulsion 
at each well U, the well depth e, and the kinetic energy 
scale of each well: 



The universal function F will behave differently in dif- 
ferent physical regimes. For example, if the wells are very 
deep, £ ^ t, the tunneling rate should be exponentially 
suppressed J ~ exp(— ae/t) with a a positive numerical 
constant. Such a case is uninteresting for the collective 
lattice dynamics of the boson system and will not be con- 
sidered here. Hence we will never write the dependence 
of the function F on t/e. We will study tunneling when 
there are a large number of bosons in the adjacent sites 
between which tunneling is occurring. This can corre- 
spond to the case where there is a very large number 
of bosons in the whole lattice or when the condensate is 
placed at two sites and its tunneling behavior monitored. 

The BEC on an optical lattice can be described by 
the Bose Hubbard Hamiltonian. Since we have a large 
number of bosons per well, there will be a BEC in each 
well with its own Bose Hubbard Hamiltonian. There will 
be a tunneling term t between adjacent wells. We then 
have 
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where V is the volume of each well ~ with a the 
interlattice spacing. We, of course, have N^V = 17. Thus 
the tunneling rate or "Josephson" current is 
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where at ar are bosonic creation and destruction op- 

k ^ 

erators with momentum k on the right well, while the 
"barred" operators correspond to the left well, e^: is the 
kinetic energy of a boson within each well. The wells 
are taken as large boxes. This approximation is valid as 
wells in optical lattices are fairly large (a ~ lOO's of nm) 
in typical lattices^]. 

As long as the bosons within each well are Bose con- 
densed, the individual well Hamiltonian can be diago- 
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nalized by the Bogolyubov methodj^]- If we keep the 
leading term in the interactions, we obtain the reduced 
Bogolyubov Hamiltonian 
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with eg = eg+2t7 Nq/Nl. The number in the condensate, 
iVo, is determined uniquely by the number of bosons in 
each well, N, and the repulsion U. We see that the rele- 
vant parameter that sets the scale for the interaction for 
the EEC is the combination UNq/Ni^. Hence the univer- 
sal dimensionless function that determines the tunneling 
current enhancement does not depend on U, Nq and 
independently, but rather 
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The reduced Bogolyubov Hamiltonian (0) is diagonal- 
ized by the Bogolyubov canonical transformation 
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where 9 is the relative phase between the superfluid order 
parameters, which is in general nonzero, and 
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happens at ~ yj 'HmU Nq / ffi N l. This defines the coher- 
ence length of the superfluid ^ ~ a/ h? /8mU Nq . 

When the system is probed at length scales L > ^ 
it responds collectively as a superfluid, while at length 
scales L < ^ it behaves like a gas of independent bosonic 
particles, thus responding dissipatively. We note that an 
inflnite system would become superfluid for any U no 
matter how small, as one can always go to a long enough 
length scale to see the system respond collectively as a su- 
perfluid. In here we do not have in mind infinite systems 
but finite wells of superfluid. In this case of finite sys- 
tems, if the kinetic energy scale of the well, T^,, is greater 
than the effective interaction, UNq/Nl, the system will 
not be a superfluid as the linear size of the system, V^^^, 
is less than the coherence length of the superfluid. In- 
dependent of whether the bosons are superfluid in each 
well, the system could still be a superfluid within the 
whole lattice. 

For the full lattice, when the kinetic energy scale is 
set by tunneling between sites (e^ ~ tk^a^), the coher- 
ence length is f ~ a^/NiJ/ANoU . The size of the lattice 
is, of course, L — N^a and the condition for superflu- 
idity (C < L) is U/t > 1/4NoNl. We see that either 
for a very large system or very large number of bosons, 
it becomes arbitrarily easy to become superfluid in the 
sense that an arbitrarily small U will order the boson 
fluid into a superfiuid at low enough temperature. Con- 
versely, as we make the system small, or we reduce the 
number of bosons, it becomes harder to order into a su- 
perfluid, which will show up experimentally as increased 
dissipation rates. 

Even though a high enough U can order the system 
into a superfluid, too high a U will make it into a Mott 
insulator|9Hl(ll. The transition to the Mott insulator can 
be continuous^ or discontinuous 1^. The continuous 
transition occurs when there is a superfluid^] in each well 
and U prevents Josephson tunneling between wells. The 
continuous transition only happens for a commensurate 
number of bosons per site. When the bosons in the wells 
are not superfluid within wells the transition is discontin- 
uous and happens irrespective of whether the number of 
bosons per site is integer or not [13. The second transi- 
tion always happens when the number of bosons per site 
or U/t are suffciently small as in these cases ^ > a. 

We now proceed to evaluate the current operator 



being the dispersion of the quasiparticle-like excitations 
of the system with e^: ~ h^k^/2m. As k ^ the dis- 
persion becomes phononic, which leads to dissipationless 
flow, i.e. superfiow, by the Landau argument We 
note that atomic BECs are superfluids as they have a 
nonzero sound speed and reduced scattering at long 
wavelengths 0, @|- The dispersion at large k becomes 
ti^k'^/2in, which leads to dissipative flow. Notice that 
the crossover between dissipative and nondissipative flow 
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where Ht — ~t (^ciio-k + ^I'^fc) tunneling part 

of the Hamiltonian and T denotes time. We note that 
the ground state of the reduced Hamiltonian |SJ) when 
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i = is an outer product of the Bogolyubov ground states 
for each weh, which we denote as our "vacuum" state 
|0). The excited states of such a system are denoted by 
\n) and are constructed by applying any number of Bo- 
golyubov creation operators to the vacuum state. The 
energies of such states, will be the sum of the en- 
ergies of each of the operators. When i ^ there will 
be tunneling between the two wells, which can be calcu- 
lated accurately from linear response theory as long as 
the number of atoms within the wells is large. For such a 
calculation the only relevant excited states are two quasi- 
particle states, one in each well, with opposite momenta, 
i.e. \n) = &t6^ -|0), with energies £'„ = 'i.E^. This cal- 
culation was first done by B. D. Josephson for Cooper 
paired fermionic superconductors We will follow his 
calculation closely. When the wells are superfluid, the 
tunneling current is dominated by coherent tunneling of 
Bogolyubov correlated pairs. 

In order to calculate in linear response theory[l^, we 
turn on Ht adiabatically from time T = — oo, start- 
ing from the ground state of the noninteracting sys- 
tem at such a time. If we expand the wavefunc- 
tion in terms of the noninteracting eigenstates \^) = 
a„(T) exp (— i£^„r/fi.)|n), we obtain the well known 
formula a„(T) = Snfi-iexp {iE„T/H){n\Ht\0)/{iEn + e), 
where the vacuum state energy Eq has been chosen to 
be 0. The parameter e was introduced to control the 
adiabatic switching of the interaction. The causal limit 
e ^ 0"*" should be always understood. Our interaction 
is turned on fully at time T = and the wavefunction 
of the interacting system is given by |^(0)) = |0) + \6(p) 
where [l^ 
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The Josephson current is given by (\I'(0)| j|5'(0)) 
(0| J|(5(^)-|- complex conjugate. With 
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analogous to the fermionic superconducting case|llj. We 
note that if the phase difference between the order pa- 
rameters, 9, is there is no Josephson tunneling. Since 
UNo/Ej:NL, we obtain for the enhancement fac- 



tor F = Jh/t 



In order to estimate this quantity we approximate 
the quasiparticle dispersion as phononic, E'i ~ 
2UNoh'^k'^ /itiNl, and we cut off the high momenta at 
fee of the order of the inverse coherence length. Convert- 
ing the sum into an integral 
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Note that this integral is infrared divergent but the di- 
vergence gets cut-off by the inverse linear well size V^^^. 
We finally obtain for the enhancement 
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In order for the system to be superfluid at all within 
a well we must have V^^^ > ^, otherwise there will 
not be Josephson tunneling per se and the calculation 
of tunneling will have to be done for noninteracting 
bosons trapped in wells. This condition implies that 
UNo/T^Nl > 1. Therefore, F > tNL/lGn^UNo. For 
typical lattices^, t is of the order of around kHz's. Nl 
is about 50, the number of bosons is about 5000, so that 
the number of Bose condensed atoms per well iVo N 
and thus Nq/Nl ^ 100. Values of U are more uncer- 
tain, but they range from tens to a hundred Hz, so that 
t/U > 10. We finally obtain F > 2. Therefore, when we 
have few hundred bosons per site, Josephson tunneling 
will be enhanced from the single particle tunneling value 
t/h. Of course, for other values of the lattice parameters, 
it need not be enhanced. In extreme cases, it can be 
enhanced quite a bit. Since the logarithm in expression 
(|17|l can never be negative due to the precondition for 
Josephson tunneling to be possible at all {V^^^ > S,), F 
is single valued. We then see that measurements of the 
Josephson "enhancement" F can be used to estimate the 
values of the on-site repulsion U if the other experimental 
parameters are known. 

We now study the case when the linear size of the 
wells is too small {V^^^ < 0- In this case the inter- 
action cannot order the system into a superfluid, and its 
only effect is a renormalization of the mass of the bosons. 
Hence the system is described by the Hamiltonian @ 
with U — and m the renormalized mass of the bosons. 
There is no Josephson tunneling per se as the system is 
not superfluid. On the other hand, there can be coherent 
quantum tunneling of the (nonsuperfluid) boson liquid 
between wells. 

The Hamiltonian (0J with [/ = is diagonalized by the 
canonical operators 
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to yield 



we obtain the tunneling current 
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The Heisenberg equations of motion lead to the following 
time dependent operators A^(T) = exp[—i{e^—t)T / h] 
and Aj^{T) = A^exp[— z(eg + <)T/?i]. The common phase 
factor —ie-^T/h will be dropped as only phase differences 
affect the physics. From these time dependent operators 
one immediately obtains the operators into which 
and their conjugates develop under Hamiltonian evolu- 
tion: 
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The "barred" and "unbarred" operators are orthogonal 
raising and lowering operators as unitary Hamiltonian 
evolution is a canonical transformation thus preserving 
the commutation relations of the operators at time zero. 

We can now calculate the coherent tunneling current 
between two wells. We take the right well to have {N + 
M)/2 bosons and the left well with (iV - M)/2 bosons. 
The initial wavefunction is then 



where the periodicity is an artifact of having a two site 
system. The proportionality of the tunneling current on 
the initial number difference of bosons between the sites, 
M, is not an artifact. Hence the tunneling current could 
be quite large if the number of bosons is very different 
between the different sites. The current can also be zero 
if the number of bosons is equal between wells. The fact 
that the current goes to zero for short times T ^ is also 
not an artifact since, for short times, the system cannot 
know if it is a closed two site system or a large system. 

In the present note we studied quantum tunneling of 
EEC's in optical lattices. For the case when bosons 
within each well are superfluid, we have Josephson tun- 
neling among the wells For typical lattice parameters the 
Josephson current will be somewhat enhanced from t/h^ 
but this need not be so. The enhancement factor will 
depend on lattice parameters. When the bosons in each 
well are not superfluid, there is coherent quantum tun- 
neling of essentially free bosons between the wells. This 
tunneling is proportional to the boson number difference 
between the wells essentially due to Bose statistics. 
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where the subscript means the smallest momentum 
state into which the bosons condensed. The wavefunc- 
tion at time T, |5'(r)), has the exact same dependence 
on the operators aj(r), aJ(T) that |^'(0)) has on a\, 
aj. In order to obtain the coherent tunneling current J 
we need to evaluate {^N{T)) = (*(r)|AA^|^'(T)) with 
AA^ = ajoo — ajao counting the excess number of bosons 
in the right well. Using 



AN = cos f — j (4(r)ao(r) - 4(T)ao(T)) 

+ism i — \ (4(r)ao(r) - 4(T)ao(T)) (22) 
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